Cosmological Density Fluctuations in Stochastic Gravity 
- Formalism and Linear Analysis - 
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We study primordial perturbations generated from quantum fluctuations of an inflaton based on the formalism 
of stochastic gravity. Integrating out the degree of freedom of the inflaton field, we analyze the time evolution 
of the correlation function of the curvature perturbation at tree level and compare it with the prediction made by 
the gauge-invariant linear perturbation theory. We find that our result coincides with that of the gauge-invariant 
perturbation theory if the e-folding from the horizon crossing time is smaller than some critical value (~ | 
slow-roll parameter |~^), which is the case for the scales of the observed cosmological structures. However, 
in the limit of the superhorizon scale, we find a discrepancy in the curvature perturbation, which suggests that 
we should include the longitudinal part of the gravitational field in the quantization of a scalar field even in 
stochastic gravity. 



Oh: 



> 

m 



o 



X 



I. INTRODUCTION 

Inflation has become the leading paradigm of the early uni- 
verse not only because it solves the theoretical difficulties of 
the standard Big-Bang scenario but also because it explains 
the origin of the almost scale-invariant primordial density per- 
turbations which have been found by the observation of the 
cosmic microwave background radiation (CMB). However, 
we may still not know some fundamental parts of the infla- 
tionary scenarios or models, mainly due to our ignorance of 
physics on the very short scale. It follows that we are largely 
interested in theoretical predictions about what we can learn 
about inflation from CMB LL J,^^. 

In order to analyze some inflation model by the observa- 
tional data of CMB, it is necessary to evaluate the primor- 
dial perturbations generated during inflation. We believe that 
quantum fluctuation of an inflaton field gives the origin of 
seeds of cosmic large-scale structures. So it may be important 
to reveal how such a quantum fluctuation in the inflationary 
phase becomes classical density perturbations. So far there 
have been a few approaches fll|2fl|2ll|2l|2l|2l|2l|2l,|4i 
Ho, 51, 52, 53]. Among them, stochasticgravity may provide 
one of the most systematic approaches IslmllTllsllgl \10l [Till . 
It can describe a transition from quantum fluctuations to clas- 
sical perturbations systematically. Stochastic gravity has been 
proposed in order to describe the behaviour of the gravita- 
tional field on the sub-Planck, scale which is affected by quan- 
tum matter fields. On this energy scale, the quantum effect of 
the gravitational field may be ignored compared with quantum 
fluctuations of matter fields. Hence the gravitational field can 
be treated as a classical one. The semi-classical approach is 
justified. 

In stochastic gravity, in order to find an effective action, we 
integrate out only the degree of matter fields by use of the 
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closed time path (CTP) 111, [H E [H ■ This is one way to 
perform coarse-graining flSll . As a result, we obtain the effec- 
tive equation of motion for gravitational field under the influ- 
ence of quantum fluctuation of matter fields, including those 
non-linear quantum effects. The obtained evolution equation 
is the Langevin type, which contains a stochastic source and a 
memory term. 

Stochastic gravity looks similar to stochastic inflation, 
which was first proposed by Starobinsky ifTol l20l I2TI I22I l23i 
24, 25, 26]. In stochastic inflation, in order to discuss the evo- 
lution of long wave modes of an inflaton scalar field, which 
play a crucial role in cosmological structure formation, the 
inflaton field is split into two parts: superhorizon modes and 
subhorizon ones. The long wave modes are affected by quan- 
tum fluctuations of the short wave modes. Then the effective 
equation of motion for the in-in expectation values of super- 
horizon modes can be derived by integrating out only subhori- 
zon modes. In most works on stochastic inflation, the scalar 
field is discussed in a homogeneous and isotropic spacetime, 
and superhorizon modes and subhorizon ones interact with 
each other through a self-interaction of the scalar field. 

In the last few years, some parts of quantum fluctuation 
of the gravitational field have been taken into account, by 
adopting a gauge-invariant variable as a canonical one to be 
quantized. This gauge-invariant variable includes not only 
a scalar field but also a longitudinal part of the gravitational 
field |27, 28, 29, 30]. Rigopoulos and Shellard have derived 
a non-linear evolution equation for long wave modes of the 
Sasaki-Mukhanov variable. The equation incorporates full 
non-linear dynamics on large scale. However, the contribution 
from quantum fluctuations of the short wave modes is evalu- 
ated based on the linear perturbation equations. Then only the 
effect from tree-level short wave modes is taken into account. 

In stochastic gravity, on the other hand, although the way 
of the course-graining is different from that in stochastic infla- 
tion, the evolution equations in stochastic gravity are similar, 
and both equations can describe the transition from quantum 
fluctuations to classical perturbations. Furthermore, stochas- 
tic gravity can incorporate nonlinear quantum effects of a 
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scalar field. This gives another advantage in stochastic grav- 
ity, as we will show below. 

When we discuss only the leading part of the curvature per- 
turbation (Q by means of linear analysis, in spite of mutual 
difference at the level of microphysics, it implies that most in- 
flationary models may be compatible with the observational 
data. This is because this gauge invariant variable C,, which 
is directly related to the temperature fluctuation of CMB, be- 
comes constant in the superhorizon region ll3ll [32I I33I1 . In 
order to make a difference between many inflationary models, 
it is necessary to subtract more information from the observ- 
able. For thi s purpo se, non-Unear effects have been studied 
intensely HSlIlllllllllllliSlillillil. Among such 
approaches to non-linear effects, stochastic gravity would be 
well-suited to compute the loop corrections induced from in- 
teraction between a scalar field and the gravitational field. It is 
because the CTP effective action includes also the non-linear 
effect of quantum fluctuations of a scalar field, and the effec- 
tive equation makes it possible to discuss the loop corrections 
to primordial perturbations. 

Roura and Verdaguer have applied the formalism of 
stochastic gravity to analyse the evolution of the primordial 
density perturbation lfT2i Il3[l . They discuss the evolution of 
Bardeen's gauge-invariant variable $ under the approxima- 
tion of de Sitter background spacetime. They showed that 
the order of the amplitude for $ is the same as the predic- 
tion by quantization of the gauge-invariant variable. How- 
ever, even at the tree level, so far, the time evolution of the 
curvature perturbation ( in the superhorizon region, which is 
important to compare with the CMB observation, has not been 
sufficiently discussed in stochastic gravity. Hence, before we 
will discuss the non-linear effects, it may be better to reformu- 
late how to evaluate the primordial perturbations in stochastic 
gravity. Since stochastic gravity includes additional effects 
such as a coarse-graining procedure to describe the transition 
from quantum fluctuations to classical perturbations, it is not 
trivial whether stochastic gravity predicts the same evolution 
of the primordial perturbations as that of the gauge-invariant 
approach. In this paper, we will consider the evolution of 
the curvature perturbation in a uniform density slicing. We 
find that stochastic gravity gives the same result as that in the 
gauge-invariant approach except for a limited case. The time 
evolution is characterized by the ratio of the Hubble horizon 
scale to the physical scale of fluctuation. We find that the am- 
plitude of ( in stochastic gravity deviates from the prediction 
of the gauge-invariant perturbation theory only when the ratio 
becomes nearly equal to zero. We discuss why two quantiza- 
tion procedures do not produce the same result. We will dis- 
cuss loop corrections both to scalar and tensor perturbations 
in |44]. 

The paper is organized as follows. In Sec. [Ill we briefly 
summarize the quantization of gauge-invariant variables. In 
Sec. Uni we present the basic idea of stochastic gravity, and 
consider the basic equations in stochastic gravity, i.e., the 
Einstein-Langevin equation, which describes the time evo- 
lution of the gravitational field affected by a quantum scalar 
field. Then we discuss the perturbations of the Einstein- 
Langevin equations around an inflationary background space- 



time. In Sec. IIVI we evaluate the correlation function of ( and 
compare the result with the prediction in the gauge-invariant 
approach. Finally, in Sec. [Vl we discuss the reason why the 
prediction in stochastic gravity does not agree with that in the 
conventional linear theory in the superhorizon region. 

Throughout in this paper, we consider a single-field infla- 
tion, which Lagrangian is given by 



(1.1) 



where 4> is an inflaton scalar field and V{4i) is its potential. 
To characterize the slow-roll inflation, we use two slow-roll 
parameters: 



(1.2) 



where H ^ a/a and = SttG are the Hubble expansion pa- 
rameter and the reduced gravitational constant, respectively, 
and V^^ = (PV/dcf)'^. As time variable, we adopt the con- 
formal time, t(< 0), and represent the time derivative by a 
prime. 



II. GAUGE-INVARIANT PERTURBATIONS 

In this section, we briefly summarize a conventional ap- 
proach for a quantization of an inflaton field, in which the 
gauge-invariant variable is used as a canonical variable. For 
more detailed explanations, please refer to the papers f35il45i 
l46ll . Throughout this paper, we follow the notation for per- 
turbed variables defined in |47]. 

In a linear perturbation theory, each momentum mode de- 
couples in the basic equations. Therefore, it is sufficient to 
consider only one mode with momentum k in the perturbed 
equations. As for the gauge condition, there are two conve- 
nient choices of time slicing to evaluate primordial perturba- 
tions. One is the slicing such that the spatial curvature per- 
turbation vanishes, 31 = 0, which is called the flat slicing. In 
this gauge choice, the equation of motion for fluctuation of a 
scalar field, is described as 



ip'} + 2J{^'j + + a''V^^)ipf 
+ [2a^V^A - + k<jg)] : 



0, 



(2.1) 



where (fif, A and dg are perturbed variables of a scalar field, 
a lapse function and a shear of a unit normal vector field to 
a constant time hypersurface, respectively, and 'K = a' /a = 
aH is a reduced Hubble parameter. Note that along with fluc- 
tuation of a scalar field, there appear metric perturbations such 
as A and (Tg|47]. In fact, the square bracket term in Eq. ( 12.11) 
represents fluctuation of the gravitational field. Using the per- 
turbed Einstein equations 



1{{A' + kcjg) + 2(J{' + 2'K'^)A 



2K^a'^V4,(p (2.2) 



'K{A' + k(jg) = -kV/ [0^0 + 5^0' (3 -e)] 



(2.3) 
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Eq. ( 12.11 ) can be rewritten by means of one variable, Lp^, as 



where 



Va, 



^(3 



0} 



(2.4) 



(2.5) 



The contribution from fluctuation of the gravitational field is 
put together into the effective mass term. If we ignore the con- 
tribution from metric perturbations, together with V^^, which 
are suppressed by the slow -roll parameter e and 7/y , Lpf obeys 
the evolution equation for a massless field in de Sitter space- 
time, i.e.. 



(2.6) 



This is why this slicing is preferred in order to consider the 
evolution of perturbations in subhorizon region. In de Sitter 
space, we easily find a two-point correlation function for a 
massless scalar field. In order to determine positive frequency 
modes, we have to impose initial conditions. When we con- 
sider only linear perturbations, it is sufficient to discuss sub- 
horizon modes on an initial time (t^), which are important for 
structure formation. In the subhorizon limit (|A:t| ^ 1), it 
is appropriate to impose that this scalar field behaves as if it 
were a free field in Minkowski spacetime, whose mode func- 
tions are given as 



1 



^ — ikr 



(2.7) 



at the initial time of inflation. Then the two-point function for 
|fcr| ^ 1 is obtained as 



(27r)3j(fc+p) 



2k^ 



(2.8) 



The other convenient gauge choice is the comoving slicing, on 
which the energy flux vanishes, — 0. It is preferred when 
we discuss the evolution of perturbations in the superhorizon 
region. Since the curvature perturbation in this slicing, 3lc, is 
related to tp j.- as 



H 

Jla = T^/ 



(2.9) 



the evolution equation for 31^ is obtained from Eq. ( |2.4| | and 
Eq. ( 12.51 ) as 



where 



z = — = sgn </>) . 

H K 



(2.11) 



Note that this is the exact equation for 5ic- To derive this equa- 
tion, neither the slow-roll approximation nor the long wave 



approximation is imposed. When we consider the large-scale 
limit of fc 0, Jlc includes a constant mode function. Using 
the approximation 



z 
z 



1 + 3e — rjv 



(2.12) 



we find the equation for 3?c as 



„ 2(l + 3e-77y) 2 



(2.13) 



If we ignore the time evolution of the slow-roll parameters, 
this equation is solved by Hankel functions. We find general 
solution as 



(2.14) 



where x = —kr and a — 3/2 + 3e — rjy. x denotes the ra- 
tio of the horizon scale to the physical size of a perturbation 
with momentum k. Since both Hankel functions (Ha ^ (x) and 

H^ ' {x)) behave as X " when x <^ 1, we find that He ap- 
proaches to constant in the superhorizon region. Hence, it is 
sufficient to evaluate a two-point function of !Kc around the 
horizon-crossing time, i.e. 



(3l,(fe,r,)3le(fe',,7')> = [^y{vfik,7j)pf{k',r,')) 



(2.15) 



where the suffix k represents that the variables are evaluated at 
the time when the fluctuation mode with k crosses a horizon. 
It follows from Eq. ( 12.15b that the amplitude of the metric per- 
turbation on given comoving scale, fc, is determined by the en- 
ergy density (or Hubble parameter H) and by deviation of the 
equation of state from de Sitter vacuum, which is described by 
the slow -roll parameter e at the time of horizon crossing. Tak- 
ing into account the fact that the Hubble parameter decreases 
and the slow-roll parameter increases as inflation goes on, the 
description Eq. ( 12.151 ) shows that inflation provides the red- 
tilted spectrum, which agrees with the CMB observation. 



III. STOCHASTIC GRAVITY 

In this section, first we shortly summarize the basic points 
on stochastic gravity. The basic equation in stochastic grav- 
ity describes the evolution of the gravitational field, whose 
source term is given by quantum matter fields. The effective 
action is obtained by integrating matter fields with the CTP 
formalism. We find the Langevin type equation, i.e, the so- 
called Einstein-Langevin equation which is analogous to the 
equation of motion for the Brownian particle. In Brownian 
motion, the deterministic trajectory is influenced and modi- 
fied by stochastic behaviour of the envkonmental source. It 
is worth noting that this Langevin type equation is well-suited 
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not only to helping us understand the properties of inflation 
and the origin of large-scale structures in the Universe, but 
also to explaining the transition from quantum fluctuations to 
classical seeds. 

The Langevin type equation is characterized by the exis- 
tence of stochastic variables and memory terms. The stochas- 
tic variables represent quantum fluctuations of the matter 
fields. Under the Gaussian approximation, its statistics is de- 
scribed by the two-point function Nabc'd' {xi,X2), which is 
called a Noise kernel. The Noise kernel is determined by the 
imaginary part of the CTP effective action. The CTP effec- 
tive action with a coarse-graining includes such an imaginary 
part. The effective equation of motion is also derived from 
this effective action. It is usually a differential-integral equa- 
tion representing the non-Markovian nature. The integral part, 
which is called the memory term, depends on the history of the 
gravitational field itself. In fact, its integrand is composed of 
fluctuations of the gravitational field 6 gab and the two-point 
function Habc'd'{xi, X2), which is called the Dissipation ker- 
nel. The Dissipation kernel is determined by the real part of 
flie CTP effective action. In Sec. ImBl and Sec. HITCl we wiU 
explicitly evaluate the Noise kernel and the Dissipation ker- 
nel. We also analyze the perturbations of the above Langevin 
type equation, i.e, the Einstein-Langevin equation. In Sec. IIVI 
using the results obtained in this section, we will evaluate the 
correlation function of the primordial perturbations. 



A. Einstein-Langevin equation 



tion: 



iUix)) = , 

{^ab{Xl)^c'd'ix2)) =Nabc'd'{xi, X2) , (3.1) 

where the bi-tensor Nabc'd' {xi , a;2) is the Noise kernel which 
represents quantum fluctuation of the energy-momentum ten- 
sor in a background spacetime, i.e., 

Nabc'd' {X1,X2) = ^Re[Fabc'd' {xi,X2)] 

= ^{{Tab{xi) - {fab{xi)),fab{x2) " {fab{x2))}}[g] , 

(3.2) 

where {X, Y} = XY + YX, g is the metric of a background 
spacetime, and the bi-tensor Fabc'd' {x, y) is defined by 

Fabc'd'{xi,X2) = {fab{xi)fc'd'{x2))[g] 

-(fab{xi))[g] (t,d'ix2))[g].0.3) 

The expectation value for the quantum scalar field is evalu- 
ated in the background spacetime g. The Noise kernel and the 
Dissipation kernel correspond to the contributions from inter- 
nal lines or loops of the Feynman diagrams, which consist of 
propagators of the scalar field and do not include propagators 
of the gravitational field as internal lines. Including the above- 
mentioned stochastic source of £_ab, the effective equation of 
motion for the gravitational field is written as 



When we consider an interacting quantum system, which 
includes the gravitational field on the sub-Planck scale, we 
may expect that quantum fluctuation of the matter fields dom- 
inates that of the gravitational field. In stochastic gravity, we 
assume that the gravitational field is not quantized but classi- 
cal because we are interested in a sub-Planckian scale. How- 
ever, it is important to take into account fluctuation of a gravi- 
tational field which is induced through interaction with quan- 
tum matter fields. In order to discuss such dynamics of the 
gravitational field, the CTP formalism is helpful. 

The effective action for the in-in expectation value of the 
gravitational field is derived by integrating matter fields. In 
I18J, Martin and Verdaguer derived the effective equation of 
motion based on the CTP functional technique applied to a 
system-environment interaction, more specifically, on the in- 
fluence functional formalism by Feynman and Vernon. This 
CTP effective action contains two specific terms, in addition 
to the ordinary Einstein-Hilbert action, describing the induced 
effects through interaction with a quantum scalar field. One 
is a memory term, by which the equation of motion depends 
on the history of the gravitational field itself. The other is a 
stochastic source £^ab, which describes quantum fluctuation of 
a scalar field. The latter is obtained from the imaginary part 
of the effective action, and as such it cannot be interpreted as 
a standard action. Indeed, there appear statistically weighted 
stochastic noises as a source for the gravitational field. Under 
the Gaussian approximation, this stochastic variable is charac- 
terized by the mean value and the two-point correlation func- 



G^''[g + 5g] = K' {T^") n[g + Sg] + 2^ 



(3.4) 



where Sg is the metric perturbation induced by quantum fluc- 
tuation of matter fields and stochastic source ^ai is character- 
ized by the average value and the two-point correlation func- 
tion Eq. ( 13.11 ). 

Note that this equation is the same as the semiclassical Ein- 
stein equation except for a source term of stochastic vari- 
ables £^ab, which represents quantum fluctuation of energy- 
momentum tensor of matter field. Furthermore, the expec- 
tation value of energy-momentum tensor includes a non-local 
effect as follows. It consists of three terms as 

{T''')R[g + Sg] = {f^'{x))[g] + {f('^^'m,Sg]{x))[g] 



d'vV^W'%]{x,y)5gM + 0{5g^) , 

(3.5) 

where t(i)''^ is the linearized energy momentum tensor and 
j^abcd jjjg Dissipation kernel. It depends both on the grav- 
itational field and on a scalar field. The evolution equation 
for a scalar field also depends on the gravitational field. As 
a result, the expectation value of the energy-momentum ten- 
sor depends directly on the spacetime geometry and indirectly 
through a scalar field. When we perturb a spacetime with met- 
ric g to that with {g + Sg), two different changes appear in 
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the r.h.s. of Eq. ( 13. St . The second term represents the di- 
rect change, which is expressed in terms of fluctuation of the 
gravitational field 5g as 

(3.6) 

where Sip"^ is defined by 

^ (VoVVoV + 7'' V.^Vj V') [g] (3.7) 

We can neglect this term safely on the sub-Planck scale be- 
cause this term is smaller by the order of (kH)^ than the 
previous term. To derive this expression, we have used the 
background evolution equation for a scalar field. 

While the third integral term in the r.h.s. of Eq. ( I3.5l l repre- 
sents the effect from the indirect change and is characterized 
by the Dissipation kernel, which is given by 

Habc'd' {X1,X2) = Hj^bl'd'i^i^^^) +^it;M'(2^1'2;2) (3.8) 
^abc'd'i^i^^^) = -^lni[Sabc'd'{xi,X2)] (3.9) 

Hitc'd'i^i^^^) = ^Im[Fahc'd'(a;i,a;2)] , (3.10) 
where Sabc'd' ixi,X2) is defined by 

Sabc'd' ixi,X2) = {T*Tab{xi)Tc'd'ix2))[g] ■ (3.11) 

T* denotes that we take time ordering before we apply the 
derivative operators in the energy momentum tensor. As 
pointed out in |8], only if the background spacetime g sat- 
isfies the semiclassical Einstein equation, the gauge invari- 
ance of the Einstein-Langevin equation is guaranteed. Hence, 
in this paper, to guarantee the gauge invariance, we assume 
the background spacetime satisfies the semiclassical Einstein 
equation. 

The Einstein-Langevin equation describes fluctuation of the 
gravitational field, which is induced by the quantum fluctua- 
tions of a scalar field. The Einstein-Langevin equation, given 
by Eq. (13.41 ) includes two different sources. One is a stochas- 
tic source £^ab, whose correlation function is given by the noise 
kernel. From the explicit form of a Noise kernel ( I3.2l i. we find 
that represents the quantum fluctuation of the energy mo- 
mentum tensor. The other is an expectation value of the en- 
ergy momentum tensor in the perturbed spacetime {g + Sg), 
which includes a memory term. The integrand of a memory 
term consists of a Dissipation kernel and fluctuation of the 
gravitational field. In order to investigate the evolution for 
fluctuation of the gravitational field, it is necessary to evalu- 
ate the quantum correction of a scalar field and calculate the 
Noise kernel and the Dissipation kernel. 



B. Perturbed Einstein-Langevin equation 

Since the CTP effective action contains the full quantum ef- 
fect of matter fields, it is expected that we can also deal with 
non-linear quantum effects such as loop corrections by means 
of stochastic gravity. Before we consider such non-linear ef- 
fects we shall discuss linear perturbations, which correspond 
to a tree-level effect in the Feynman diagram. Note that in a 
linear perturbation theory in the homogeneous and isotropic 
universe, three types of perturbations, scalar, vector, and ten- 
sor perturbations are decoupled from each other Then, we 
discuss each perturbation independently. 

We are interested in density perturbations as seeds for struc- 
ture formation. Then, it is sufficient to consider only scalar 
perturbations. In this paper, we choose the metric as follows: 

ds^ = -a^{T){l + 2AY)dT^ -2a^{T)^^YjdTdx^ 

ji 

+a^{T)-f,^dx'dx\ (3.12) 

where a and jij are a scale factor and the 3 metric of max- 
imally symmetric space, and A and (k/J-C)^ are the lapse 
function and the shift vector, respectively. The scalar pertur- 
bations can be expanded by a complete set of harmonic func- 
tion Y^) on the three-dimensional space and Yj is defined by 
Yj = —k^^Yy. In this gauge choice, the spatial curvature 
perturbation vanishes. We will discuss the evolution of the 
following variable: 

C^Ys^f^ (3-13) 

where Sf = Sp/p is density perturbation. This variable C is 
gauge-invariant, and turns out to be a curvature perturbation in 
a uniform density slicing. In the classical perturbation theory, 
the energy conservation law implies that this variable is con- 
served in the superhorizon region for a single-field inflation 
ll3ll [32i I33I1 . ( is directly related to a gravitational potential at 
the late stage of the universe and in turn to the observed CMB 
fluctuations. 

The density perturbation in the present slicing is given by 

STl EE -pSfY 

= Sgoc{f''{x))[g]+goc{{f('^°n<l>[9lSg]{x)) 

-2 J d^y^^)H^^^'^[g]ix,y)Sgdeiy) + 2ei ■ 

(3.14) 

Since the background energy-momentum tensor is given by 

= g^\f\{x))[g] = -pg"^ , (3.15) 

the direct contribution from a change of the gravitational field 
is described as 

{f^'^°"m,5g]{x)) = + '- + 0{{KHf)}pAY. 

(3.16) 
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With these two relations ( |3.14| l and ( |3.16t . the density pertur- 
bation in flat slicing is written as 



2s 2 

6 p 



(3.17) 



where we have defined the density perturbations of the 
stochastic source and of the memory term as follows: 



5p„ 



d'^xe' 



ikx 



X Hoac'd' [g] (x, y)9'''^' g'^'^'Sge' f (y) 



-.900 e'{x) 



(3.18) 



Since the contribution from the memory term, 6pm , contains 
the past history of metric perturbations, it seems very difficult 
to deal with such a term. However, as we will show in Sec. 
IIV Al we can evaluate the memory term, rewriting the integral 
equation into a differential equation. 

The Hamiltonian constraint equation gives a relation be- 
tween the gauge-invariant variable A and the density pertur- 



bation S f as 



Using it, we eliminate A in Eq. ( |3.17| i, and find 



1 / fc \2 



(l ~l)^f = -piSPi + Spm) + 0(( W) . 



(3.19) 



(3.20) 



Hence, in the superhorizon region, the two-point function for 
Sf is expressed in terms of four correlation functions of Sp^ 
and 6p,n, i.e, 



+ i^P,nki^)^Pipi^)) + {SPik('^)^P"^pi^)) 



+ {^P,nk('^)^P"^pi'^)) 

Here we have used the relation 

y(r) = (l-|)p + 0(p i>,Hf) 



(3.21) 



(3.22) 



The four correlation functions in Eq. ( I3.21l i are described by 
the noise and the Dissipation kernels. We present the details 
in AppendixlAl The correlation function for Sp^ is given from 
the Noise kernel as follows: 



('^Pefc(^i)%p('r2)) 



K^eViV2S{k+p)Re G+(ri, ra) 



(3.23) 



where G\{ti, T2) is the Wightman function in momentum 
space, which will be given by the mode function. We have 
denoted X{Tj) as Xj, where X{t) is some function of a con- 
formal time T. 



Since the Noise kernel is not a time-ordered correlation 
function, it should be described by the Wightman function in 
momentum space G\ (ri , r2). In order to derive this relation, 
we have to assume that 



d 
a Oti 



Re 



(t>__d_ 
a dT2 



■Rc 



« 



a^^^nVRe G+(ti, T2) 



(3.24) 



where a^^^ = V,f,/{KV). The Wightman function with the 
initial condition imposed in Sec lIIICI satisfies this relation in 
slow-roll inflationary universe. 

While the density perturbation of the memory term is ex- 
pressed as follows : 



(3.25) 



Sp,„,ki^) ^ -2eK'V{T) / dri(ai) 



As mentioned in Appendix |A] since contributions from sub- 
horizon region will oscillate, we can ignore them in the inte- 
gral of Eq. ( 13.25b . Hence the lower bound of the time inte- 
gral in Eq. (13.25b may be given by the horizon crossing time, 
Tk ~ — 1/fc. Since the memory term depends on the past his- 
tory, the integrand of the memory term contains 6pm itself. 
In Sec. IIVI rewriting this form to a differential equation, we 
evaluate contributions from the memory term. 

Then, once the Wightman function, G\{t, T2), is deter- 
mined, we can evaluate the correlation function for 5 f and ( 
from Eq. ( 13.21b . (13.23b . and (13.25b . In the next subsection, im- 
posing the appropriate initial condition, we will find the mode 
functions, and in turn the corresponding Green function. 



C. Wightman function 

Here we consider the Wightman function in momentum 
space. 



G+(ri, r2)^V^fc(ri)V^(T2), 



(3.26) 



where i/;/^ (r) is the mode function of a quantum scalar field 
in the inflationary universe. It satisfies the wave equation 



(3.27) 



We solve this equation under the slow-roll condition. Intro- 
ducing a new variable as ipj^ir) = a(r) iPi^{t), this equation 
is rewritten as 

^fc(T) + [e~{2~e- 77v(3 - e)}D{^]i^f,iT) = 0, (3.28) 
where we have used the relation 



In an inflationary stage on the sub-Planck scale, we can ne- 
glect the term whose magnitude is smaller by the order of 
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{kHY than that of the leading term. We also ignore non- 
linear terms with respect to the slow-roll parameters. So we 
do not take into account time evolution of the slow-roll param- 
eters. Under these assumptions, the equation for -ip becomes 



dx 



1 



2 + 3{e-rjv) 



Tp{x) = 0, (3.30) 



where wehaveusedJ{~ — 1/[(1 — £)t]. The general solution 
for this equation is given by the Hankel functions as 



a;2 



13 



(3.31) 



where 0^ = 9/4 + 3(£ — ?7y), with two arbitrary constants C 
and D. It implies 



X2 



CH};\x)+DHf\x) 



(3.32) 



We can assume that the mode functions should have the same 
form as in Minkowski spacetime, i.e., 



1 



2k 



(3.33) 



when the wavelength is much shorter than the horizon scale, 
i.e., at very early times of the universe. This fact may be true 
in the present gauge rather than the comoving gauge. Then 
the mode function and the Wightman function in momentum 
space are given by 




Hf\x) 



(3.34) 



fjfcln, T-2j- Hp (Xi) (X2), 

4 aia2 

(3.35) 

where x = —kr. Setting a.; — 1, the scale factor a(T) is given 

by a(r) = ijr^ ' Using this expression, the Wightman 
function is rewritten as 

_ Tryri T2 / Tir2 \ i+g 



Gtiru r2) = 



H^;\x^)Hf{x2) 



fc3 V rf J 
Here we have used the relation 



(3.36) 



(3.37) 



In order to compute the correlation functions for 6p^ and Spm, 
it is sufficient to consider the evolution of the Wightman func- 
tion in the superhorizon region. The behaviour of G^t (ti , T2) 
in the superhorizon region is given in Appendix iB] Substitut- 
ing the approximated expression ( |B5t for the real part of the 



Wightman function into Eq. (|3.23t , we find that the correla- 
tion for the density perturbation of stochastic variable 6p(^ is 
expressed as 

.Pf i(.+rt.V„,4|^<£ii|)!l^ ,338) 

4 fc'^ xf^ 

Similarly, substituting the approximated expression ( IB6I ) for 
the imaginary part of the Wightman function into Eq. ( 13.251) . 
we find that the density perturbation for the memory term is 



2 /■! 



3 afc2 

/3 / a; \/5 



dxi{ai)^Xi 



X / \a;i 

2(l+e) ,1 



fc2 



dxiXi 



-3e 



IV. 



{(?) -(f) +^''..*('.)). 

(3.39) 



CORRELATION FUNCTIONS IN STOCHASTIC 
GRAVITY 



As shown in Eq. ( |3.21t . the correlation function for the 
density perturbation in flat slicing consists of four correlation 
functions : 

T2) = (<5p^fc(Ti)%p(r2)) (4.1) 

/Cm(Tl, T2) = {Sp^k('^l)Spmp{T2)) (4.2) 
/t<(ti, T2) = ('5p„fc(Tl)%p(T2)) (4.3) 
/mm(Tl, T2) = {Sp^f.{Ti)Sp.mp{T2)) . (4.4) 

Here we have not explicitly shown the momentum depen- 
dence in those expressions just for simplicity, /t ^ has already 
been given in Eq. ( 13.38b . i.e., 

T2) ~ ^-^ b{k +p) eViV2- 



k^ X 



2e 



(4.5) 



The other three functions contain the contribution from the 
memory term. In the next subsection, IIV Al we present these 
three correlation functions. Then, these functions determine 
the correlation functions for 5f (or Q, which also gives the 
correlation function for the curvature perturbation in uniform 
density slicing. This gauge- invariant variable C is related to 
the fluctuation of the temperature of CMB. 



A. Memory term 

The memory term represents the non-Markovian nature of 
the effective equation of motion. The evolution equation de- 
pends on the past history, showing its non-local nature. The 
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memory term describes interesting phenomena, such as dissi- 
pation. It plays an important role in non-equilibrium systems. 
Especially, in stochastic gravity, the memory term represents 
the indirect dependence on the gravitational field in the energy 
momentum tensor This point has already been emphasised in 
Sec. IIII A I The analysis, however, bothers us due to the com- 
plexity of the integral equation. Here, rewriting this integral 
equation into the differential equation, we evaluate the contri- 
bution from the memory term. 

The time evolution of the density perturbation for the mem- 
ory term 5pm is described by Eq. ( |3.39t . It determines the 
correlation function for Sp^ in the integral form. The cor- 
relation function between 5p^ and 5pm satisfies the integral 
equation: 



/r<(a;, X2) = eH"^ — 



3+e 2(l+£) 



dxiXi 



-3e 



(-) 



(4.6) 

Taking derivatives of this equation twice, we find a differential 
equation as 

d 



_1 r -2-2(£-w)_l 

dx I dx 



frnd^,X2) I 



2l3ex 



-4-4e+r(i 



We can simplify it to the second order differential equation 
with a source term, i.e., 

■^fmdx,X2) - -^frnd^'^^) 

+ ^^^T^frad^^ ^2) = X2). (4.8) 

Since Eq. ( |4.6l l is the integral equation, we have to impose the 
initial conditions, i.e.. 



fmdl,X2) = —fr„dx,X2) 



0. 



x = l 



(4.9) 



The time x — 1 corresponds to the horizon crossing time. 
Using two independent solutions for the homogeneous equa- 
tions, which are the power-law functions 



3-(e+r)i/) 



f2{x) 



(4.10) 



we find the solution for Eq. ( |4.8l l with initial conditions ( |4.9l l 



fmiix,X2) = e 



^ dxi 

Xl 



(-) -(-) (4.11) 



Substituting /tj given in Eq. (14.5b into Eq. ( 14. lib , we obtain 
the correlation function for 5pm and Sp^ as 

(27r)3 



/mc(-ri,r2) 



6{k+p) eViV2^^{x,X2Y^ 



fc3 

{x^'-l) + \{{x,f-l} 



(4.12) 



where we have used the fact that the Hubble parameter at the 
initial time is related to that at the horizon crossing time as 
Hfxl'^^ ~ Hi Note that /^^(ri, ra) = /^^(ra, n). 

As for the correlation function for Spm, we have the integral 
equation. 



fmmix, X2) = SH 



3 , J 2(l+e) 



dxix^ 



-3e 



isobtained 



( — ) - ( — ) ■ iUmixi, X2) + fmmixi, X2)] ■ 

(4.13) 



It is converted into the equivalent differential equation, 

, ^ Jmm\X^X2j , Jnim\X^X2j 

dx'^ X dx 



3(£ + rjv) 



fmm (*^; -^2 ) 



%Umix,X2) 



X 



with the initial conditions, 

d 



/mm(l,a;2) 



dx 



fmm {x^ X2 ) 



= 0. 



(4.14) 



(4.15) 



Using the time evolution of f^m{xi, X2), we determine the 
time evolution of fmm{xi, X2) as follows: 



[hd^,X2)+fmdx,X2)]. (4.7) fmm{xi,X2) 



1 dx 



3-(e+r)v) 



X f^mix, X2) 



S(^^{k+p) eV,V2^{x,x2r- 



fc3 

(^r-i) + |{(^i)'-i} 

(^2"^-l) + |{(^2)^-l} 



(4.16) 



Here, we find the following specific property of the memory 
term. The leading term of the equal time correlation function, 

fmm 

(t, t), is given by 



fmm {^1 



X X 



2j7v 



(2^)3 T.2K'^fc 



{x-^-l) 



(4.17) 



After this fluctuation crosses the horizon scale, until 
— elogx « 1, x~'^ is nearly equal to unity. Hence the contri- 
butions from the memory term are suppressed by the slow-roll 
parameter e, compared to /{^(t, r). In other words, as long as 
the e-foldings from the horizon crossing time to time r, which 
is given by Nk~tT — — log x, is smaller than 1/e, the contri- 
bution from the memory term is negligible in the correlation 
function for S / (or Q. 

However, once the e-foldings Nk^r gets larger than 1 /e, 
the term x^'^ becomes much larger, and eventually the con- 
tribution from the memory term dominates the direct contri- 
bution from the stochastic variable /^{(t, t). In this case, the 
correlation function is approximated as 



fmm{^: 



^^5{k+p)sV'^x'^^^-^\ 
4 P 



(4.18) 
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B. Correlation functions 



V. SUMMARY AND DISCUSSIONS 



Substituting the correlation function for 5p^ and 5pm, given 
by Eq. (|431l, Eq. (|4TT2] | and Eq. (|4TT6] | into Eq. dO), we 
determine the correlation function for the density perturbation 
in flat-slicing, 5f, as 



-P) £ 
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(4.19) 



The density perturbation in flat-slicing, Sf, is related to the 
curvature perturbation in uniform density slicing, (, by Eq. 
( 13.131 1. We also find the correlation function of ( as 



(Cfc(T)Cp(T)) 



(27r)3 K^Hi 



(27r)= 



^-2 IJ2 



(4.20) 



Note that after the horizon crossing time, until the e-folding 
Nk^T approaches to l/jryy — sl^ the correlation function for 
C is approximated as 



(Cfc(r)Cp(r)) 



(27r) 



-S{k +p) 



-'Hi 



(4.21) 



Then, in this region, the amplitude of the curvature perturba- 
tion C is constant. However, once the e-folding Nk^r exceeds 
l/|?7y — e|, it will deviate from the constant value. 

In order to compare our result with the prediction from 
quantization of the gauge-invariant variable summarized in 
Sec. |II] it is helpful to show the relation 



B = 



9(1 



A: \2 



(4.22) 



which is derived from the momentum constraint in the uni- 
form density slicing. This relation implies that in the super- 
horizon region, the uniform density slicing agrees with the 
comoving slicing||48l]. Hence, the correlation function for the 
curvature perturbation C in the gauge invariant perturbation 
theory is given by that for 'Re, i-C-, Eq- ( 12.151 ). Comparing 
Eq. ( 14.201 1 to Eq. ( |2.15t , we find that as long as the e-folding 
from the horizon crossing time is smaller than 1 / 1 ?7v — e | , our 
result is the same result as that in the gauge-invariant pertur- 
bation theory. However, when the e-folding Nk^r becomes 
larger than 1 / 1 t^v/ — £ | , the predictions by these two approaches 
are inconsistent. 

In the next section, we discuss the reason why this deviation 
appears, and why the amplitude in stochastic gravity does not 
predict a constant evolution if the curvature perturbation is in 
the the superhorizon region. 



In this paper, we have considered the evolution of the pri- 
mordial perturbations in a slow-roll inflationary universe. Pri- 
mordial density perturbations, which play a crucial role at the 
later structure formation stage, are generated inside the hori- 
zon from the quantum fluctuation of a scalar field in the be- 
ginning of inflation, and then stretched out to the superhorizon 
scale. 

However, since the observable quantity through the CMB 
observation is classical, in order to constrain an inflation 
model by the observation, it is necessary to consider a tran- 
sition from quantum fluctuations to classical perturbations. 
Coarse-graining is the well-suited way to understand such a 
transition. Hence we have discussed cosmological perturba- 
tions in the inflationary stage in stochastic gravity. Based on 
the naive expectation that quantum fluctuation of the gravi- 
tational field is sufficiently small at the sub-Planck scale, we 
assume in stochastic gravity that contributions from the Feyn- 
man diagrams in which the gravitational field propagates as 
an internal line are negligibly small. Then, when we compute 
the CTP effective action, we integrate out only the dynam- 
ical degree of freedom of a scalar field. As a result, from 
this coarse-grained effective action, the evolution equation for 
the gravitational field is derived, which is called the Einstein- 
Langevin equation. 

Although this effective equation makes it possible to dis- 
cuss non-linear quantum effects such as the loop corrections 
of a scalar field, in this paper, we have first considered the 
tree-level effect, which is the leading contributor to the pri- 
mordial perturbations. We have derived two-point correlation 
functions for density perturbations and the curvature pertur- 
bations. We find that our results are consistent with those in 
the gauge-invariant perturbation theory unless the e-folding 
from the horizon crossing time exceeds some critical value 

(i/\w~£\y 

However, as pointed out in the previous section, the curva- 
ture perturbation ( does not keep constant in the superhori- 
zon region if the e-folding gets beyond the critical value, in 
contrast to the prediction of the gauge-invariant perturbation 
theory. We expect that such a deviation appears due to neglect 
of quantum fluctuation of the gravitational field. In stochastic 
gravity, we do not integrate out the degree of freedom of quan- 
tum fluctuation of the gravitational field in the path integral of 
the CTP effective action. In this sense, quantum fluctuation 
of the gravitational field is neglected, although we include the 
fluctuation of the gravitational field induced by the quantum 
fluctuations of matter fields. It seems that this missing part 
induces the deviation in the superhorizon region. 

In order to understand this point more clearly, it is helpful 
to reconsider the time evolution of the mode functions, tpj^, 
which equation is given by Eq. ( |3.27| i. We emphasize that 
this equation represents the wave equation for quantum scalar 
field in the background inflationary spacetime. 

Comparing Eq. ( I3.27l i with Eq. ( 12. Il l for fluctuation of a 
scalar field in flat slicing {(p /), we find that the equation for Lp f 
contains the fluctuation of the gravitational field. Such a con- 
tribution plays a crucial role in keeping the curvature perturba- 
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tion constant in the superhorizon region. In a gauge-invariant 
theory, the curvature perturbation in comoving sHcing, 3?c is 
related to (yS/ as 



H 

'Xc = T-(p/, 



(5.1) 



and the equation for 3^c is given by Eq. ( 12.1 11 1. 

In comoving slicing, 3lc is the only one dynamical degree of 
freedom, and we have a single equation for it without neglect 
of any perturbed variables. The equation guarantees that the 
curvature perturbation is constant at a superhorizon scale. 
It means that only when all contributions of perturbations in- 
cluding fluctuation of the gravitational field are taken into ac- 
count, C=const in the superhorizon region is guaranteed. If 
any part of the contributions is ignored, it implies the devia- 
tion from the constant evolution of C, in the superhorizon re- 
gion. 

In fact, when we transform il) into ifjc = {H/ (t>)ip, imitating 
the relation between 5(f} / and 'Re, the equation for -ipc is written 



- 3 - e - 



^^fc(r)=0 (5.2) 



Unlike the case of 'Jlc, because of the existence of the term 
proportional to e in the coefficient of ijj^ ^, V'^ y[j.=constant is 
no longer the solution for the superhorizon scale {x <C 1), 
although the deviation is very small in a slow-roll inflation. 

Finally, we consider the reason why neglection of quan- 
tum fluctuation of the gravitational field has influenced the 
behaviour of perturbations in the superhorizon region. In this 
discussion, first we have to distinguish fluctuation of the lon- 
gitudinal mode of the gravitational field induced by matter 
fields from fluctuation of gravitons. 

The latter may be very small and can be ignored at least 
on the sub-Planckian energy scale, while the former may not 
be so, because it is difficult to distinguish fluctuation of a 
scalar field from that of the longitudinal part of the gravita- 
tional field. In fact, if we change the gauge condition, they are 
mixed up. Hence when we quantize a scalar field, it may be 
natural to include the longitudinal gravitational field as well 



H54I1 . In the present calculation based on stochastic gravity, 
we ignore both quantum fluctuations of the gravitational field. 
For the sub-Planck scale inflation, however, even in stochastic 
gravity, we should include contributions from fluctuation of 
the longitudinal mode which couples to the matter field, as we 
discussed above. We may have to improve our formulation 
for density perturbations in stochastic gravity. However, we 
should emphasize that the deviation in the present approach 
from the gauge-invariant approach becomes large only when 
the e-foldings from the horizon crossing time to the definite 
time exceed l/|e — r}v\- For example, if we assume that 
\e ~ rjv\ is around 0.01, then it implies that this deviation 
appears only when the e-foldings from the horizon crossing 
time to the end of inflation Nk^^ exceeds about one hundred. 
Taking into account that e-folding Nk^e for today's Hubble 
horizon scale is about fifty, we may conclude that neglect of 
the longitudinal modes will influence only on a larger scale 
than the observed region today. 

In stochastic gravity, as we mentioned, we can calculate the 
loop corrections, because the CTP effective action includes 
also the non-linear effect of quantum fluctuations of a scalar 
field, and the effective equation makes it possible to discuss 
the loop corrections to primordial perturbations. We will dis- 
cuss it in iHl. 
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APPENDIX A: NOISE KERNEL AND DISSIPATION KERNEL 



In AppendixlAl we present the Noise kernel and the Dissipation kernel, which are given by the bi-tensors, Fabdd' {xi , X2) and 
Sabc'd' {xi , 2:2 )■ These expressions are given from integration of the CTP effective action in terms of fluctuation of a scalar field. 
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As seen in ||44J . the lowest-order part of these bi- tensors can be given by 

- Vabiaif KVxa^^\2^2 (<5°, G + 6°^, G+ , + r^c'd' G+ ,) 

- Vc'd'{a2f KV2a^i\i^i {5\ G+ + 6\ G% + lUb G+ ) 
+ ^:6€','(«ia2)'«'V^i^^2a«a«G+ 



(Al) 



Sabc'd'{xi,X2) = aia20102 { 5''j\,{lGF);bd' + 5''jd'{lGF);bc' + 5\5l,{lGF),ad' + 5\5%{lGF);ac' 

+ 'nc'd'iS"ai'i-GF):bO' + S°i,iiGF):aQ') + VabiS^c' i^G f) -Qd' + <5 (^^F ) ;0c' ) + VabVc' d' [iG f) -00'} 

- Vabiaif KVia[^'' a2<i>2 (-5"^' {iG f) -d' + S°d' {iGF),c' + Vcd' {iGF),o') 

- Vc'd'{a2Y KV2a^2^ai4>i (5°^ {iGF)-b + S\ {iGF)-a + Vab (*Gf);o) 



+ V^bVl'd'i^i^^f k'ViV2 a['^a^i\iGF) 



(A2) 



where a*^^^ = V^/kV ~ sgn(V0)V2£- There appear the two different Green functions such as G"*" = G+(xi, X2) and Gf — 
Gf{xi, X2)- The former is the Wightman function, and the latter is the Feynman function. To compute the correlation function 

Sp^ and Sp,n, we have to know Fqoo'o' {xi, X2), Faoo'i'{xi,X2), Sqqq'o' (xi, X2), and Sooo'i' {xi, X2), which are written by 



F°oo'o'i^i^ ^2) = -{a2fOr,Or2G+{xi, X2) 



00'; 



(xi, a;2) = -a2(/>20TiG+(xi, a;2):(' 



'5'"oo'0'(a;i> 2^2) = -(02)^ 6'(ti - r2)0^i0^2G+(a;i, a;2) + 0{t2 - Ti)Or^Or2G+ {x2,Xi) 

S°0Q'1'{XI, X2) = -a24>2 0{ti - T2)OtiG^{xi,X2);1> + 0{t2 - Ti)0r^G^(x2,Xi).i> 



where we have defined the operator Or as 



(j) d 



a(r) dr 

Using Eq. (IA3I )-( IA6I) . we shall evaluate the noise and Dissipation kernels in order 



(A3) 
(A4) 

(A5) 
(A6) 



(A7) 



1. Noise kernel 

Substituting the expression Eq. ( |A3t into Eq. ( |3.2| l. we obtain 

{iV^i)fii'{x2)) = nY^,{x^,x2) = iRe[^^oOo' (^1,2:2)] = io.,0.,Re[G+(xi, ^2)]. 
Hence we find the correlation function for 5p^ as 

{5p^k{r^)5p^p{T2)) - j d^x, J d^X2e-^^-^^ e-^P-^He,{xi)e'A^2))'^^'^ 



4 

(2^)3 



(5(fc+p)0,,0,,Re^G+(Ti, T2) 
5{k+p)eK^ViV2Re\G\{Ti, T2) 



(A8) 



(A9) 



In the last equality, we have approximated the operator as Or — a^^^ nV, because the time derivative of the Wightman function 
drG^ is proportional to the slow-roll parameter rjv (see AppendixlEt. 



12 



2. Dissipation kernel 

Substituting the expression Eq. iA3i and Eq. (|A5l l into Eq. (13.101 1 and Eq. ( 13.91 ), we obtain 



H 



(2)0 



,o,(a;i, xa) = H^^^^ ^Q,g,{xi, X2) + H)t'" ^^,^,{xi, X2) 



(2)0 



00 



1 



-Im 



.(2)0 



00 



,Q,{xi, 2:2) + S'^^^"oo,o,(xi, X2) 



(«2)^ 



^^(Tl -T2)0.,0.Jm G+(a;i,a;2) 



Similarly, substituting the expression Eq. (IA4b and Eq. ( IA6b into Eq. (13.101 1 and Eq. (13.91 ), we find 



H 



(2)0 



X2) = a;2) + i?^"^ (0:1, X2) 



(2)0 



-Im 



0.2 



F 



(2)0 



YiAxi-, X2) + S'^^^"oo,;,(a;i, a;2) 



= --^ 4'2 9{ti- T2)Or^lni G^{xi,X2)-l' 

In order to evaluate 5p_^ ^ (t), it is convenient to introduce the Dissipation kernel in momentum space as 

Habc'd'{Tl, T2, p, q) = S'^^^ip + q) Habc'd'in, T2, 

= J d^x, J d^X2 e-^P-^^e-^'i ''^Habc'd'{x,, X2). 
Especially, the components, (0, 0, 0', 0') and (0, 0, 0', V), are given by 



H 



(2)0 



OO'O 



,o,(ri, T2, k) ^ 



{02? 



(2^)30(ti - T2) Or.OrM G\{ti, T2) 



H^'^\o'i'i^i^ r2, k) = <P2 {2nfe{n - T2)^ O.Jm G+(ri, T2) 



Using this expression, 5p^ (r) is written as 



(ai)^[^"oo^d'(^, n, k'')Ai.{Ti) + i^.. 00'/' V, jr^ 



(2), 



dTi {aiY 



G\iT, n) 



(MO) 



(All) 



(A12) 



(A13) 



(A14) 



where we have neglected the contribution from tensor perturbations which is expected to be much smaller than that from scalar 
perturbations and 'I'^), and for the second equality, we have used the relation ( |3.19t . 

We can neglect the first term in Eq. ( IA14l i as follows. In superhorizon region (xi < 1), this term is suppressed by the factor 
of (/c/!Ki)^. Around horizon crossing time {xi ~ 1), this term is suppressed by the slow -roll parameter Neglecting this term 
and using the approximation of Or — a^^^'nV, we find the density perturbation from the memory term Spm as 



rfri (ai)'*(5^j^(Ti) OrOrJm G\{t, ti) 



2eK'V / dn (ainSp^ki^i) + <5p,„fc(ri)} Im G+(r, n) 



(A15) 



The lower bound of this integration is given by the horizon crossing time Tk = —1/k because in the subhorizon region (a;i > 1), 
the Hankel function appearing in G\{t, ti) oscillates, and then no accumulative contribution exists. 
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APPENDIX B: WIGHTMAN FUNCTION IN THE SUPERHORIZON REGION 

In this section, we will consider the evolution of the Wightman function in the superhorizon region, which is given by 



G^in, T2) = -—— Hy{xi) Hy{x2). 

Since the Hankel function is expressed in terms of the Bessel function as follows : 



^ Vf" sin/37r' 

The Bessel function is expanded by a power series of x in the region of < a: < 1, i.e.. 



i~irix/2) 



2m 



to! r(/3 + TO + 1) 

Then we obtain the asymptotic behaviour of the Hankel function in the superhorizon region as 



Re[HP{x^)Hy>{x2)] 



(2) 



lm[H,^'\x,)H,^'\x2) 



sin/37rr(l - /3) 
cos(/3- 3/2)7r 



1 



sin2/3^ r(i - /3)r(l + /3) 



Using these formulas, the asymptotic behaviour of the Wightman function in the superhorizon region is given by 

:/TlT2V 



Re[G+(ri, r^) 



Mr ^/TiT2 



Im[G+(Ti, T2)] ~ M 



(1 — e)^ aia2 

^/TlT2 



0102 



{X1X2)-^ 
P /.TiN/3 



MMix^X2)U^y {X^X2)-'' 



\XiJ Va;2/ 



where the coefficients Mj. and A/j are defined by 

M 



22/3 ^ 

[sin/?.r(l-/3)P = 2+^('' '^^^ 



TT cos(/3-3/2)7r 
4 



sin^ r(l - /3)r(l + /?) 6 



0(e, r)v)- 



(Bl) 



(B2) 



(B3) 



(B4) 

(B5) 
(B6) 

(B7) 
(B8) 



In order to evaluate the correlation function of 5p^, it is also necessary to compute the time derivatives of the real part of the 
Wightman function, which are given by 



dr,Re[G+{Ti, T2) 



G+(ti, T2), 

Tl 

dr.,Re[G+in, r2)]^^G+(ri, r2) , 

T2 



Vv 
T1T2 



(B9) 
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